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Abstract 

We show how a particular spatial structure with a buffer globally stabilizes the chemostat dynamics 
with non-monotonic response function, while this is not possible with single, serial or parallel chemostats 
of the same total volume and input flow. We give a characterization of the set of such configurations that 
satisfy this property, as well as the configuration that ensures the best nutrient conversion. Furthermore, 
we characterize the minimal buffer volume to be added to a single chemostat for obtaining the global 
stability. These results are illustrated with the Haldane kinetic function. 
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1 Introduction 

The chemostat has been introduced in the fifties as an experimental device to study of microbial growth on 
a limiting resource [29[ 132) . It is also often used as a mean to reproduce situations where (limiting) nutrients 
are fed to micro-organisms, typically in a liquid medium, such as in natural ecosystems [161 15] or anthropized 
environments |23j . More generally, the chemostat is largely used as a scientific investigation tool in microbial 
ecology [2TII45]. 

The mathematical model of the chemostat has been extensively studied (see e.g. [3D]) and used as a 
reference model in microbiology [33] , microbial ecology |10) or biotechnological industries such as the waste- 
water treatment industry [7J. However, in many applications, the assumption of perfectly mixed chemostats 
has to be relaxed. In the eighties, the gradostat, as an experimental device composed of a set of chemostats of 
identical volume interconnected in series, has been proposed to represent spatial gradient [25] . For instance, 
it has been used to reproduce marine environment [18] or to model rhizosphere [11] . and has motivated 
several mathematical studies (43] [20l [9] [37J [47] EH [T7j [4TJ [12] . Similarly, an interest for series of bioreactors 
appeared in biochemical industry, with tanks of possibly different volumes to be optimized [26j [El H] [15j IB] . 
Although island models have been proposed in ecology since the late sixties [27J , relatively few studies have 
considered non-serial interconnections of chemostats [38] . In natural reservoirs such as in undergrounds or 
ground- waters, a spatial structure with interconnections between several volumes is often considered, each 
of them being approximated as perfectly mixed. Those interconnections can be parallel, series or built up in 
more complex networks. To our knowledge, the influence of the topology of a network of chemostats on the 
overall dynamics has been sparsely investigated in the literature However, the simple consideration of two 
different habitats can lead to non-intuitive behaviors [321 1301 GS3 [22] and influence significantly the overall 
performances [31] [13] . 
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It is also well-known since the seventies that microbial growth can be inhibited by large concentrations of 
nutrient. Such inhibition can be modeled by non-monotonic response functions [l]|3] and lead to instability 
in the chemostat [U HH1 HI] • Several control strategies of the input flow have been proposed in the literature 
to globally stabilize such systems [5] HH GUI 133] but the ability of a spatial structure to passively stabilize 
such dynamics has not been yet studied (in [35] a general structure of networks of chemostats is considered 
but with monotonic growth rates, while in |44) non-monotonic functions are considered but for the serial 
gradostat only). 

The present work considers non-monotonic response functions with a particular interconnection of two 
chemostats of different volumes, one being a buffer tank. To our knowledge, this spatial structure, that 
is neither serial nor parallel, has not yet been considered in the literature. The idea is to decouple the 
residence time of microorganisms in two vessels such that the wash-out equilibrium is repulsive in both 
tanks. We prove that this is possible with such a configuration, while any serial, parallel or single tank 
structures with the same total volume exhibits bi-stability. This result brings new insights in microbial 
ecology for the understanding of the role of spatial patterns in the stability of bio-conversion processes in 
natural environments, where natural buffers can occur, such as in soil ecosystems. It has also potential 
impact on the design of robust industrial bio-processes. 

The paper is organized as follows. Section [5] presents the hypotheses and the buffered configuration, 
comparing with serial and parallel interconnections. Section [3J studies the multiplicity of equilibriums and 
provides a complete characterization of the set C of configurations that have a unique positive equilibrium. It 
is then shown that this equilibrium is globally asymptotically stable. Section 2] analyzes the performance of 
the buffered chemostat, characterizing first the best configurations among C in terms of nutrient conversion 
at steady state, and then the minimal volume of a buffer to be added to a single chemostat to make the 
interconnected system globally stable. Finally, in Section [5] numerical simulations illustrate the results on 
the Haldane function. 



2 General considerations 

We consider the chemostat model with a single strain growing on a single limiting nutrient. The system 
is fed with nutrient of concentration Sin with flow rate Q. The total volume V is assumed to be constant 
(i.e. input and output flow rates are supposed to be identical). When the concentrations of nutrient (or 
substrate) and biomass, denoted respectively S and X, are homogeneous, as it is the case in perfectly mixed 
tanks, the system can be modeled by the well-known equations: 

X = v(S)X-^X, 

where /i(-) is the uptake function and Y the yield coefficient of the transformation of nutrient into biomass. 
Changing the units in which time t, growth /x(-) and biomass X are measured, one can assume without any 
loss of generality that D — Q/V — 1 an Y = 1 (replacing t by t/D, fi(-) by Dji(-) and X by YX). 

In this work, we consider specifically strain growths that present an inhibition, described by the following 
assumption. 

Assumption Al. The function /i(-) is C°°([0, +oo)) and such that /i(0) = 0, fi(S) > for any S > 0. 
Moreover there exists a number S > such that \i is increasing on (0, S) and decreasing on (S, +oo). 



For instance, the Haldane function [T] 



= K + S + P/K, (2) 
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fulfills Assumption Al (see Figure [TJ. Classically, we consider the set 

A = {S> | n(S) > 1} (3) 

that plays an important role in the determination of the equilibriums of the system. Under Assumption Ai, 
the set A is either empty of an open interval that we shall denote 

A=(A_,A+) , 

where A+ can be equal to +oo. 

We recall from the theory of the chemostat model [3^ the three kinds of solutions flow of the dynamics 
(fTj) under Assumption Al, depending on the parameter Si n . 

Proposition 1. Assume that Hypothesis Al is fulfilled. 

- Case 1: A = or A + > Si n . The wash-out equilibrium Eq — (Si n ,0) is the unique non negative 
equilibrium of system (Qp. Furthermore it is globally attracting. 

- Case 2: Si n > A+. The system has three non-negative equilibriums E- = (\-,Si n — A_), E + = 
(A+,iSi n — A+) and Eq = (Si„,0). Only £L and Eq are attracting, and the dynamics is bistable. 

- Case 3: Si n G A. The system (QP has two non negative equilibriums £L = (A_,5 jrl — A_) and 
Eq = (Si n ,0). E- is globally attracting on the positive quadrant. 

In case 2, the issue of the growth can change radically depending on the initial condition. Throughout 
the paper, we shall consider uptake functions /i(-) and values of Si n that fulfill the conditions of Case 2, so 
that the system exhibits bi-stability in the classical chemostat model ([1}. 

Assumption A2. A_ < A + < Si„. 




Figure 1: Graph of the Haldane function and illustration of Assumption A2. 

The question we investigate in this paper is related to the assumption that the vessel is perfectly mixed, 
and to the role that a spatial structure could have on the stability of the dynamics. We shall consider spatial 
configurations with the same input flow and residence time than the perfectly mixed case, that is to say with 
the same total volume V and input flow Q. Under Assumption A2, one can check that having a volume 
V split in several smaller volumes Vi interconnected in series or in parallel, assuming that each of them 
is perfectly mixed, leads necessarily to a global dynamics having bi-stability or the wash-out as the only 
equilibrium in at least one of the tanks. 

- In the series connection, the dynamics of the first tank of volume V\ is given by equations ([lj where 
V is replaced by V\ < V. Its dilution rate is then equal to Q/Vi, that is greater than Q/V = 1. 
According to Proposition [T] only Cases 1 or 2 can occur. 
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- In the parallel connection, the dynamics of each tank of volume V, and flow rate Qi is given by equations 
(P) where V and Q are replaced by Vi and Qi. Denote r*j = Vi/V and on — Qi/Q, and note that one 
has Si r « = Si a i = Then, the dilution rate Di in the tank i is equal to oii/ri. According to 
Proposition [TJ a necessary condition for having a single attracting equilibrium in each tank is to have 
Di < 1 for any i, that contradicts J2i r i = Si a i = 1- 

In the present work, we study a particular spatial configuration with a asymmetry created by two inter- 
connected volumes, one of them serving as a buffer (see Figure[5|), that we shall call the "buffered chemostat" , 
to be compared with the "single chemostat" . V\ and Vi are respectively the volumes of the main tank and 




Ql + Q 2 

Figure 2: The buffered chemostat. 



the buffer, and Q\ and Q2 denote the input flow rates of each tank. We assume that each vessel is perfectly 
mixed. Straightforwardly, the dynamical equations of the buffered chemostat are 



Si 

s 2 
x 2 



QiSi n + Q2S2 — QS\ 
Vi 

Q2X2 — QX\ 



-V(S 2 )X2 + 



QiSin — Q2S2 



(4) 



fi(S 2 )X 2 



Q2X2 

Vn 



v 2 



Note that the limiting case V\ = consists in a by-pass of the volume V2 with a flow Qi 



In the following sections, we shall consider configurations with the same total volume V = V\ + V2 and 
input flow Q — Qi + Q2, to be compared with the single chemostat (Vi = V and V2 = 0). In Section|4l we 
shall also consider configurations with a fixed volume Vi = V and study the benefit of adding a buffer of 
volume V 2 , under a constant total input flow Q = Qi + Q2- 



3 Study of multiplicity of equilibriums 

Given fi(-) and Si n such that Assumptions Al and A2 are fulfilled, we describe the set of all possible 
configurations such that Q = Qi + Q2 and V = VI + V2 (with Q/V = 1) by two parameters r € (0, 1) and 
a > defined as follows 

Vi Q2 

T = — ry = 

V (l-r)Q 
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Dynamics (U) can then be written in the following way 

a(l - r)(S 2 - 5i) + (1 - a(l - r))(S in - Si) 



5i = -//(5i)Xi 



r 



Ai = ^(^ijAiH , ^ 

•52 = -[J.(S 2 )X 2 + a(S in - S 2 ) , 
X 2 = v{S 2 )X 2 - aX 2 . 

One can easily see that an equilibrium (S^, X*, S 2 , X^ ) of dynamics (jSJ is the solution of the following 
equations: 

1 + ^ (l - «f^fl) = ^(^) or S i* = S ™ > ( 6 ) 

A'* = iSj n — S* , (7) 
/i(S*) = a or S* = S in , (8) 

X 2 = Si n — S 2 . (9) 

Due to the cascade structure of the model (Q} , the study of the dynamics of the second reactor can be done 
independently of the first one. Depending of the value of a, the three cases depicted in Proposition [T] for the 
single chemostat are possible. 

Under Assumptions Al and A2, one can straightforwardly check from equations (JSJ) and © that there 
exists an unique positive equilibrium (S 2 , X$) in the second tank exactly when a belongs to the set (0, fJ-(Si n )] . 
For any number a £ (0,//(Sj n )], we then denote S 2 (a) = S 2 £ (0, Si n ) the unique solution of the equation 

M (S 2 *(a))=a. (10) 

We shall consider the family of hyperbola H a r that are the graphs of the functions 

a (\ 1 , 1 ~ r {■, S in -S${a)\ 

9W(s) = H 1-a — (11) 



parametrized by a € (0, /Lt(5j n )] and r € (0, 1). From equations © and ([7]), a positive equilibrium (S*,X*) 
of ([5]) has to fulfill precisely 

that is to claim that S* is the abscissa of an intersection of the graph of /i(-) with the hyperbola H a r . To 
each solution S± corresponds an unique X$ = Si n — S*. We define the family of sets 

11(a) = {re (0, 1) | 3ls £ (0, S m ) s.t. <t> a , r (s) = fi(s)} (12) 

parametrized by a £ (0,//(Sj n )]. The set C of pairs (a,r) such that dynamics ([5]) admits an unique positive 
equilibrium is then defined by 

C = {(a,r)\ae(0,fi(S in )],ren(a)} . (13) 

For convenience, we shall consider the set of s at which the hyperbola H a ,r is tangent to the graph of 
the function /i(-) and is locally on one side (that amounts to have as a local extremum of the function 

<?W(-) - /«(■) at s ) : 

f f d n (b a r d n U 1 1 

S a r = < s £ (X-,Si n ) s.t. min < n £ IN I — r~ i ~( s ) —, — ( s ) r is even and larger than 1 > (14) 

ds n ds n 
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and define the number 

5(a) — aS^ia) + (1 — a)Si n . (15) 
We consider two subsets of values of r such that the hyperbola H a r is tangent to the graph of /i(-). 

R~(a) = {re (0, 1) | 3s e S a , r with (s - 5(a)) (A + - 5(a)) < 0} , (16) 

R + (a) = {re (0, 1) | 3s e S a , r with (s - A+)(A+ - 5(a)) > 0} . (17) 

We state now our main result about the multiplicity of equilibriums of system ([5]) and give a charac- 
terization of the sets IZ(a) defined in (|12l) . depending on the sets R~(a), R + (a) and their interlacing. In 
Section [H this result is applied to the Haldane function ([2]). 



Proposition 2. Assume that Hypotheses Al and A2 are fulfilled. For any a £ (0,/i(5i„)] and r £ (0, 1) the 
dynamics {5p admits a positive equilibrium with 5* such that 

(5(a) - 5*)(A+ - 5(a)) > . (18) 

The set i? + (a) is non empty, and the set R~ (a) is not reduced to a singleton when it is not empty. One has 

(0, min R + (a)) when R" (a) = , 



K{a) = 
Furthermore, 



(0,mini?+(a)) H (0,1) \ [min R~ (a) , max R~ (a)] when R~ (a) ^ . ^ 19 ^ 



- for any r £ (min i? + (a), 1) ; there exists at least two equilibriums such that (5(a) — 5*)(A+ — 5(a)) > 0, 
and at least four for r in a subset of (mini? + (a), 1) when i? + (a) is not reduced to a singleton, 

- when R~(a) is non empty, for any r 6 (mini?~(a),maxi?~(a)), there exists at least three equilibriums 
such that (5(a) - 5£)(A+ - 5(a)) < 0. 

Proof. Fix a £ (0,^(5i„)] and simply denote by 5£ and 5 the values of 5£(a) and 5(a). For each 
r £ (0, 1), we define the function 

fr(s) = <f>a,r(s) ~ A*( s ) ■ 

A positive equilibrium for the first tank has then to satisfy f r (S*) = 0. 

One can easily check that <fi a ,r(S_) = 1 whatever the value of r £ (0, 1). 4> a . r (') being decreasing, one has 
4>a,r{s) > 1 for s < 5 and <f> a ,r{s) < 1 for s > 5. For convenience, we shall also consider the function 

^ = c 5 f ( V TTT ^ 2 °) 

S - S in + [bin ~ S)fl(s) 

that is defined on the set of s £ (0, 5 !n ) such that (5i„ — s)/i(s) ^ Sm — 5. On this set, one can easily check 
that the following equivalence is fulfilled 

f r (a) = <=► 7 (a) = r . 

From (|2T)|) . one can also write 

(Ms) - 1)^ 

7(s) = 



t> a>r (s)-l)j^:-l + fJ,(s) 
and deduce the property 

i(s) = <==> <j>' P (*)(/*(*) - 1) = Wa,r(«) - l)M'(s) ■ (21) 
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Recursively, one obtains for every integer n 



{£Z(.) = 0,p = l. {^( s )(„( s ) - 1) = (^, r ( s ) - 1)^( S ) , p = 1 ■ • -n} . 
Consequently, the set S a ^ r defined in (|14l) can be characterized as 

S a , r = |s G (A_, Si n ) s -t- j(s) = r and min |n £ IN* | ^-^(s) 7^ Oj is even| 

or equivalently 

$a,r = {s G (A_ , s.t. 7 (s) = r is a local externum } . (22) 
We distinguish now several case depending on the position of S_ with respect to A + . 



Case 1: S < A + . 



If S_ < A_, one has > and f r (S) < for any S € A. / r (-) being decreasing on [0, A_], one deduces 

that there exists exactly one solution S 1 * of f r (S) = on the interval [0, A+], whatever is r. Furthermore, this 
solution has to belong to [S_, A_]. The functions 4> r (-) and /i(-) being respectively decreasing and increasing 
on this interval, one has necessarily 7 / (<S'j[ r ) 7^ and then R~(a) = 0. 

If S > A_, one has f r (S) > for any S G [0, A_], and / r (5) < for any 5 G [S, A+]. On the interval 
I = (X-,S), the function is well defined and 7(1) = (0, 1) with 7 (A_) = 1 and 7(5) = 0. If R~{a) is 
empty, then 7(-) is decreasing on /, and for any r G (0, 1) there exits an unique 5* G I such that j(S*) = r. 
If R~(a) is not empty, property (|22[) implies that 7 admits local extrema. Let r~, be respectively the 
smallest local minimum and the largest local maximum of 7 (-) on the interval /. Note that r~ > and 

< 1 because 7(1) = (0, 1), and that one has r~ = mini?~(a) < = maxi?~(a). By the Mean Value 
Theorem, there exists exactly one solution 5 1 * of j(s) = r on the interval [0, A+] for any r ^ [ r mi r ~M\i an( i 
there are at least three solutions for r G (r^.r^). 

Consider now the interval J — (A + , Sj n ) where the function 7 (-) is well defined and positive with 7(A + ) = 
1 and lim s _ i .5 iji 7(5) = 1. We define 

r + — min{7(s) | s G J} 

that belongs to (0, 1). Then r + belongs to R + (a), and for any r < r + there is no solution of j(s) = r on 
J. Thus r + is the minimal element of R + (a) . By the Mean Value Theorem there are at least two solutions 
of 7 (s) = r on J when r > ? ,+ . When i? + (a) is not reduced to a singleton, the function 7 has at least on 
local maximum tm and one local minimum r m , in addition to r + . By the Mean Value Theorem, there are 
at least four solutions of 7 (s) = r on J for r G (r m , rjw). 

Finally, we have shown that the set R + (a) is non empty, and that the uniqueness of the solution of 
7(5'*) = r occurs exactly for values of r that do not belong to the set [mini? - (a), max i? - (a)] U [min R + (a), 1]. 



Case 2: S= X+. 



One has f r (S_) — for any r, so there exists a positive equilibrium with S* = S_. f r (S) > for any 
S e [0,A_] and the function 7(-) is well defined on I U J = (A_,5) U {S_,S in ) with 7(1 U J) = (0,1), 
7 (A_) = 1 and lim s _ ) .s , <n 7( s ) = 1- Using the L'Hopital's rule, we show that the function 7 (-) can be 
continuously extended at S_: 

y .S,. -T 1 

S ™ TiSJ -™ + (Sin ~ s)/i'(s) 1 " (Si„ - ' 



7 



Note that n'(S_) < so that 7(5) belongs to (0, 1), and we pose 

f = min{7(s) | s £ (A_,5; n )} . 

Then, for r < f, there is no solution of y(s) = r on (A_,5j„), and 5 is the only solution of f r (s) = 
on (0, Si n ). On the contrary, for r > f, there are at least two solutions of j(s) — r on (A_,5j„) and the 
dynamics has at least two positive equilibriums. 

Similarly, the function 7(-) is C 1 on (A_,5m) because it is differentiable at 5: 

(5 m -5K(5)-2^(5) 
7l - J [l-(5 in -5>'(5)] 2 

(and recursively as many time differentiable as the function //(•) is, minus one). Then f is the minimal 
element of the set R + (a), and the set R~(a) is empty by definition. As previously, when R + (a) is not 
reduced to a singleton, j(s) = r has at least four solutions for r in a subset of (r + , 1). 

Case 3: 5 > A + . 

We proceed similarly as in case 1. Note first that there exists no solution of f r (s) — on the intervals 
(0, A_) and (A+,5) whatever is r. 

On the set A, j(-) is well defined with 7(A) C (0, 1), 7(A_) = 1 and 7(A+) — 1 and we pose 

r + = min{7(s) | s £ A} 

that belongs to (0, 1). One has necessarily r + = mini? + (a), and there is no solution of 7(5*) = r exactly 
when r < r + . For r > r + , there exists at least two solutions by the Mean Value Theorem, and four for a 
subset of (r + , 1) when R + (a) is not reduced to a singleton. 

On the interval J = (5, 5j„), the function j(-) is well defined with 7 (J) = (0, 1), 7(5) = and 7(5j„) = 1. 
There exists at least one solution of f r (s) = on this interval. If R~{a) — 0, 7(-) is increasing and there 
exists an unique solution of 7(5*) = r on J whatever is r. Otherwise, mini? - (a) and maxi? _ (a) are the 
smallest local minimum and largest local maximum of the function 7 on the interval J, respectively. Then, 
uniqueness of 5* on J is achieved exactly for r that does not belong to [mini? _ (a),maxi?~(a)], and for 
r £ (mini? _ (a), maxi?~(a)), there are at least three solutions by the Mean Value Theorem. □ 

For each a £ (0,/i(5j„)], we define the number 

f(a) = supTZ(a) . (23) 
The following remark deals with the continuity of the map f(-). 

Remark 1. According to Proposition^ for any a and r such that r £ 71(a), one can define uniquely a 
number 5*(a,r) £ (0, Si n ) such that 

Q , r (5*(a,r))= M (5*(a,r)) . 
The map (ce,r) 1— > 5*(a, r) is clearly continuous and one can then consider the limiting map: 

5*(a) = lim 5*(a,r) . 

Furthermore, accordingly to Proposition^ one has 5* (a) < A + (resp. 5*(a) > A+j when 5(a) < A+ (resp. 
5(a) > \+). Consider, if it exists, a value of a, denoted by a, that is such that 5(a) = A+. Although one 
has 4> a ,r(^+) = MA+) for any r, there is no reason to have 

lim 5* (a) = A + or lim 5* (a) = A + . 

a<a,a— ¥a a>a,a— >a 

Consequently, the map a 1— > f (a) might be discontinuous at such point a. 
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We consider now pairs (a,r) 6 C and denote by S*(a,r) the corresponding value of Si at the unique 
positive steady state, that we denote E*(a, r). 

Let us consider first the (62,^2) sub-system, that is the the single chemostat model. Several results in 
the literature have shown the global stability of the positive equilibrium for this model, provided that the 
condition fi(Si n ) > a is satisfied (see for instance proofs for the non-monotonic case in [5J [H]). We 
present here a little extension of this result that shows global stability in the limiting case a = fJ,(Si n ), that 
is not covered in the statement of Proposition [1] 

Lemma 1. For any configuration (a,r) £ C and non-negative initial condition with X 2 (0) > 0, the solution 
62(4) and X 2 (t) of |5|) is non negative for any t > and one has 

lim (S 2 (t),X 2 (t)) = (S^(a),S m - S£(a)) . 

t— >+oo 

Proof. From equations ([5]) one can write the properties 

s 2 = ==► £ 2 > , 
x 2 = x 2 = , 

and deduces that the variables S 2 (t) and X 2 (t) remain non negative for any positive time. Considering the 
variable Z 2 = S 2 + X 2 — whose dynamics is Z 2 = — aZ 2 , we conclude that S 2 (t) are X 2 (t) are bounded 
and satisfy 

lim S 2 {t) + X 2 {t) = S in . 

t— f+oo 

The dynamics of the variable S 2 can thus be written as an non autonomous scalar equation: 

S 2 = (a- A*(S , 2 )(S'm - S 2 ) - ii(S 2 )Z 2 (t) 

that is asymptotically autonomous. The study of his asymptotic dynamics is straightforward: any trajectory 
that converges forwardly to the domain [0, Si n ] has to converge to S 2 (a) or to Si n . Then, the application of 
Theorem E] (see Appendix) allows to conclude that forward trajectories of the (5*2,^2) sub-system converge 
asymptotically either to the positive steady state (S 2 (a) 7 Si n — S 2 (a)) or to the "wash-out" equilibrium 
(Si n ,0). We show now that for any initial condition such that ^(0) > 0, the forward trajectory cannot 
converge to the wash-out equilibrium. From equations one can write 

(n(S 2 (r))-a)dT 



X 2 {t) = X 2 (0)eJo 
If X 2 (.) tends to 0, then one should have 

+00 

0i(5 2 (r))-a)dr = -oo (24) 

for any finite positive T. Using Taylor-Lagrange Theorem, there exists a continuous function 6(.) in (0, 1) 
such that 

M(S 2 (r)) - /i(5<„) + fJ,'(S 2 (r))(S 2 (T) - S in ) + KS m ) - a with S 2 (t) = S m + 9(r)(S m - S 2 (r)) . 
One can then write 

+00 r-\-oo r-\-oo r-\-oo 

(fi(S 2 (T)) - a)dr = / (fi(S in ) - a)dr - ^ {S 2 {r))X 2 {r)dr + n> ' {S 2 { T ))Z 2 { T )dT 
Jt Jt jt 

r + OO 1 r+cc 

> -/ fl'(S 2 (T))X 2 (T)dT- - / n'(S 3 (T))Z 2 (T)dT . 



9 



Note that 5*2 (.) tends to Si n when X2{-) tends to 0. So there exists T > such that ^(r) > S for any 
t > T, and accordingly to Assumption Al, there exist positive numbers a, b such that — //(S^t)) 6 [a, 6] 
for any r > T. The following inequality is obtained 

r+oo r+oo i 

/ (m(S 2 (t)) - a)dr > a / X 2 (r) - -|Z 3 (T)| 
leading to a contradiction with (pM|) . □ 
Then the global stability of the positive equilibrium of dynamics ([5]) can be proved. 



Proposition 3. For any configuration (a, r) 6 C, any trajectory of the dynamics f5|] with -X^(O) > 
converges exponentially to the steady state E*(a,r) in forward time. 

Proof. Let us consider the vector 

r, _ X\ + S% — Si n 

[ X2 + S2 — Si„ 

whose dynamics is linear: 




A 

The matrix A is clearly Hurwitz and consequently Z converges exponentially towards in forward time. 
Furthermore, variables S2 and X2 being non negative (see Lemma [T|), one has also from <j5j) the following 
properties 

Si = ==> Si > , 
Xi = => Xi > , 

and deduces that variables Si and Xi stay also non negative in forward time. The definition of Z allows us 
to conclude that variables Si, Xi, S2, X2 are bounded. 

From equations ([S]) , the dynamics of the variable Si can be written as an non-autonomous scalar equation: 

Si = (-KSi) + 1 ~ a(1 ~ r) ) (S in - Si) + a{1 - r) (S a (t) - Si) - ^(Si)Zi(t) 
\ r J r 

or equivalently: 

Si = (K, r (Si) - ^(Si))(S m - Si) + a — (S 2 (t) - S 2 *(a)) - ^(Si)Zi(t) . (25) 

r 

By Lemma [1] we know that S 2 (£) converges towards S%(a). So the dynamics (|25p is asymptotically au- 
tonomous with the limiting equation 

Si = (4> a ASi) - MSi))(S m - Si) . (26) 

Proposition [2] guarantees that S*(a,r) is the only solution of 4> a , r {Si) — y>(Si) on the domain (0,Sj n ). So 
any trajectory of ([2^)1 that converges forwardly to the domain [0, Sj„] has to converges to S*(a, r) or to Sj„. 
Along with Theorem [6] (see Appendix), we conclude that forward trajectories of the (Si,Ai) sub-system 
converge asymptotically either to the positive steady state {S*(a, r), Si n — S*(a, r)) or to (S;„, 0). We show 
that this last case is not possible. From equations ([5]), one has 

Xi = 0=^Xi = ^^l X 2 
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but from Lemma [TJ we know that X%{t) converges to a positive value and consequently X± cannot converges 
towards 0. 

Finally, we write the Jacobian matrix J* of dynamics J3]) at steady state E*(a,r) in (Z, Si, S%) coordi- 
nates: 



J = 



A 





-fi(St) 


r 


- M (S 2 *) 


-Li'(SZ)(S in - S|) 



V 

Remind the following facts: 

i. A is Hurwitz, 

ii. S% is less than S, so one has ^'(S^) < (cf Assumption Al) 

hi. Si is the only zero of <p a , r (Si) = n{Si) on (0, Si n ), the graph of (/> QjT . is not tangent to the graph of 
^ (Proposition ^ and a , r (O) > = 0, so one has <j>' ar {S*) - ^{S{) < 

One can then conclude that J* is Hurwitz. □ 



4 Study of performance of the buffered chemostat 

We first aim at characterizing among all the configurations in the set C the ones that provide the best con- 
version of the nutrient at steady state, that is the smallest value of S* (a, r). 

For convenience, we consider the function 

1>(8)=n(s)(Sin-s) (27) 

and define the number 

V>* = max ip(s) (28) 

sG[0,s] 

where s is defined by 

S = lim 5£(a) . (29) 

The number s is such that fx(s) = fJ,(Si n ) with s < Si n . 

Assumptions Al and A2 provide the uniqueness of s* realizing the maximum in (1281) , and one can then 
define the number 

a* = fjt( 8 *) . (30) 

Lemma 2. Assume that Hypotheses Al and A2 are fulfilled. For any a € (0, (i(Si n )], one has the following 
property 

- ifS(a) < X+. 

inf S*(a, r) — lira S*(a,r) , 

r^TZ(a) r— >f(a) 

- if S_(a) = A+, S*(a,r) is equal to A + whatever is r £ TL{a), 

- ifS{a) > \+, 

inf St(a,r) = lim St (a, r) = 5(a) . 
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where S_(-) is defined in \15)) . 

Proof. Fix a 6 (0, fi(Si n ))- One can check from expression (fTTj) that the map 



4> a ,r(s) is 



decreasing for s € [0,5(a)) 
increasing for s £ (5(a)), 5j., 



Along with property (fT8|) of Proposition [2] one deduces the following properties. 

i. When 5(a) < A+ and r £ 71(a), the unique positive S*(a,r) solution of (j) a , r {Si) = /i(5i) belongs to 
[0,5(a)] and r H> S*(a,r) is decreasing. 

ii. When 5(a) = A+ and r £ 71(a), S± (a, r) = A+ is the only positive solution of 4> a ,r(Si) = fi(Si). 

iii. When 5(a) > A + and r £ 71(a), the unique positive S*(a,r) solution of Q , r (5i) = m(5i) belongs to 
[5(a), Si n ] and r S*(a,r) is increasing. 

The conclusion comes then straightforwardly. □ 

Lemma 3. Assume that Hypotheses Al and A 2 are fulfilled. The following property is then satisfied. 

inf St far) = inf St(a*,r) . 

(a,r)eC relZ(a) 



Proof. Remark first from (fTOj) and (p~5|) that one has 

5(a) = 5 m - tp(S%(a)) 

and 

ae(0,^,„))^S 2 *(«) G (0,s) . 

We consider now three cases depending on ip* and A + . 

If xj}* < Sin — A+, then for any a £ (0, fj,(Si n )], one has 5(a) > A+ and according to Lemma [2] one has 

inf Si(a,r) = 5(a) . 

Then one can write 

inf 5i*(a, r) — inf 5(a) = S(a*) — inf 5i*(a* r) . 

If ip* — Si ri ~X+ then 5(a) > A + for any a £ (0, /i(5i„)], and according to Lemma[2]one has S* (a, r) > A + 
for any r £ 71(a). Furthermore, for a — a* , one has S_(a*) — A+ for any r £ 71(a). This implies the equality 

inf S*(a, r) = A + = 5(a*) = inf St(a\r) . 

(a,r)eC rER(o) 

If ?/>* > 5i„ — A + , then according to Proposition [5] there exist values of a £ (0, (j,(Si n )] such that 
S*(a,r) < 5(a) for any r £ 71(a). Then one can write the following inequality 

, fG *, u 1 i 1 ~ r f S(a) - S^(a,r) 

1-r f S( a *)-Sj(a,r) \ 
> 1 + S m -Sl(a,r) J =«W(* («.'))• 

Remind that one has (j> air (S*(a.,r)) — fi(S*(a,r)). This implies that the root 5*(a*,r) of the function 
s i — y 4> a *,r(s) — n(s) is necessarily such that 5*(a*,r) < S*(a,r), for any r £ 71(a). Finally, we obtain the 
equality 

inf St(a,r) = inf St(a k ,r). 

(a,r)eC r£K(a) 
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□ 



Lemmas [5] and [3] give the following characterization of the best configurations. 

Proposition 4. Assume that Hypotheses Al and A2 are fulfilled. The best stable configuration consists in 
choosing a — a* (or a arbitrarily close to fJ.(Si n ) if a* = (J,(Si n )) and 

- having a by-pass of the volume V with a flow rate equal to (1 — a)Q, when ij)* < Si n — A+. The output 
concentration at steady state is then equal (or arbitrarily close) to Si„ — ij)*. 

- choosing any value of r 6 lZ(a), when ip* = Si n — A + . The output concentration at steady state is then 
equal (or arbitrarily close) to A + . 

- taking r smaller and arbitrarily close to f(a), when if)* > Si n — A + . The output at steady state is then 
arbitrary close to the infimum of S* on S (that is necessarily less than \+). 



Under Assumptions Al and A2, we study now the benefit of adding to a single chemostat of volume V a 
buffer of volume V2 under a constant total input flow Q = Qi + Q2, and characterize the minimal value of 
V2/V to obtain a global stability of the positive equilibrium. Similarly to Section [3l we describe the set of 
configurations by two non-negative parameters: 

Q2 R V 2 

a = ^' p= v> 

but here one has V\ = V whatever is the volume V2. For any number a € (0, fi(Si n )], there exists an unique 
S2 (a) 6 (0, s) defined by (fT0|l and (|29j) . and consequently there exists an unique positive equilibrium in the 
second tank. The parameter a being fixed, one can straightforwardly check on equations ((U) that a positive 
equilibrium in the first tank fulfills 

<p(St) = aP(S in - S$(a)) (31) 



where the function if is defined as 



¥>(*) = (S in -s)(l-ii(s)) . (32) 



Consequently, we are looking for the smallest value of (3 such that there exists an unique positive solution 
of (|3Tj) on the interval (0, Si n ). 



Proposition 5. Under Assumptions Al and A2, there exists a buffered configuration with an additional 
tank of volume V% that possesses a unique globally exponentially stable positive equilibrium from any initial 
condition with £2(0) > 0, exactly when V2 fulfills the condition 

max (p(s) 

where tp* is defined in (28\) , with any a £ (0, fJ,(Si n )] such that 

max <p(s) < af3(Sin ~ (a)) < Si n ■ 

Proof. Let us examine first some properties of the function if on the interval (0, Si n ): 
. ip is negative exactly on the interval A, 
. tp' is negative on (0, A_) with ^(0) = S m and ^(A_) = 0, 

. = tp(Si n ) — and (p reaches its maximum m + on the sub- interval (A+, Si n ), that is strictly less 

than S in = <p(0), 
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Figure 3: Illustration of the graph of the function tp 



from which we deduce that there exists a unique solution of <p(s) — c on the whole interval (0,5; n ) exactly 
when c € (m + , 5j„) (see Figure (|3]) as an illustration). 

Consequently, the configurations (a,/3) for which there exists a unique 5*(a,/3) € (0, 5j„) solution of the 
equation (j3"Tj) are exactly those that fulfill the condition 



a(S in - S£(a)) 



</?< 



5; 



a(5j„ - 5£(a)) 



with a £ (0,/x(jSi„)]. Then, Proposition [3] with r = 1/(1 +/3) guarantees that the unique positive equilibrium 
(5* (a, /3),Si n — 5* (a, /3), 5^ (a), 5, n — 5£ (a)) is globally exponentially stable on the domain ]R^ x 1R+ x M + . 

Among all such configurations, the infimum of /3 can be approached arbitrarily close when a is maximizing 
the function 

a i y a(S ln - 53(a)) 

on (0, fx(Si n )], that exactly amounts to maximize the function ip defined in (|27[) on the interval (0, s], □ 

Remark 2. Section /ias shown the benefit of the buffered chemostat in terms of global stability of the 
system, but with a price to pay in performances when one imposes to have the same residence time (i.e. 
the nutrient concentration at steady state is larger than A_J. When adding a buffer, this is no longer true 
(i.e. the steady state necessarily exhibits a better performance than A_J: there always exists a solution 
S* 6 (0, A_) of L'Jl]) . that is unique under conditions of Proposition^ because f(0) = 5 jrl , <^(A_) = and 
a/3(5 m -5*) e (0,5 m ). 

5 Illustration and discussion 



We illustrate the results of the former sections on an non-monotonic uptake function given by the Haldane 
expression ([2]). One can easily check that for this function the set A defined in ([3]) is non empty exactly 
when the condition 

m> 1 + 2W— 

V K-I 

is fulfilled. Then, A_ , A + are given by the following expressions: 



A± = 



K^ft - 1) ± y/K-fifl-lf-AKKj 



Remind that Assumption A2 is fulfilled for values of 5j„ larger than A + . 
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Lemma 4. Assume that /x(-) is an Haldane function and that Assumptions Al and A2 are fulfilled. For any 
a £ (0, fx(Sin)], the following properties are satisfied. 

- the set R + (a) defined in |_?7| ) is a singleton, 

- for any r £ (0, 1), the set S r . a defined in is either empty or a singleton, 

- if the set R~(a) defined in \lb}) is non empty, then one has maxi? _ (a) < R + {a). 
Proof. In the case of the Haldane function, the equality <p a , r (s) = /i(s) can be rewritten as 

{S m -s-a(l- r)(S m - S$(a))(K + S + S* 2 /A7) = rfis(S in - s) . 

So S* is the root of a polynomial P of degree three, and there exists at most three solutions of (j> a , r {s) — fi(s). 
We then deduce from Proposition [2] that R + (a) is a singleton. 

Requiring to have 4> a ,r{s) = n(s) and <f>' a r {s) = fi'(s) simultaneously implies that s is solution of P = 
and P'(s) = i.e. that s is a double root of P. P being of degree three, there is a most one such solution. 
So the set S rya possesses at most one element, and this implies R~(a) D R + (a) = 0. 

When R~ (a) is non empty, we know from Proposition [5] that for r £ (mini? _ (a),maxi? _ (o;)), <f) a ,r( s ) = 
/j,(s) has at least three solutions on an interval /, and for r £ (mini? + (a), 1) at least two on another interval 
J, where I and J are disjoint. Consequently, one should have maxi? _ (a) < mini? + (a), otherwise there 
would exists at least 5 solutions of (j) a . r (s) — /i(s) on (0, Si n ). □ 

Lemma 2] implies that for any a £ (0,/i(Si n )], the number f(a) defined in (|23p is the single element of 
R + (a). It can then be determined numerically as the unique minimizer of the function 

F a (r, s) = (fi(s) - ^Asjf + (n'(s) ~ ^M) 2 
on (0, 1) x {s £ (A~, Sin) s -t- (s — A+)(A+ — S.(a)) > 0} that is, for the Haldane function: 

F ( )-{ V s 1 1 - r S ln - SZ(a) \ 2 f ^K-s 2 /Kj) 1 - r S in - S$(a) x 1 

a[r,S> {K + s + syK! r +a r S m -s J \(K + s + s 2 / 'Kjf r (S m - s) 2 

where S% defined in ([TU)) is given by the formula 



_ Kjjfi -a)- yj Kj(fl - af - Aa 2 KKi 
b2[a) ~ Ya ' 

and Sip) is defined in ([15)) . 

One can also easily check that for the Haldane growth, the function ip defined in (|27|) is increasing up to 
ip* and decreasing. Its maximum on the interval [0, Si n ] is achieved for the value 



_ yjK 2 + KS ln (l + Sin/Kij - K 

1 + S m /Kj 

Consequently, one has 

s* = min(s*, s) , 
that allows to determine the optimal value a* = £t(s*). 

The parameters given in Table [1] have been chosen for the numerical simulations. 

Figure 3] illustrates the family of functions 4> a .r(') an d the tangent property with the graph of an Haldane 
function. On Figure [SJ the domain C defined in (|13l) is drawn for different values of Si n . According to 
Remark [T] one can see that the map a <-> f (a) is discontinuous at a = a, where a is such that s(a) = A + 
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A* 


if 




A_ 


A+ 


12 


1 


0.8 


~ 0.103 


~ 0.777 



Table 1: Parameters of the Haldane function and the corresponding values of A_ , A+. 




Figure 4: Family of functions 4> a .r{ ) when 5(a) < A + (in the left) and 5(a) > A+ (in the right). 




Figure 5: Domain C of stable configurations for different values of Si 



(when it exists). On Figure |6]one can see that the two limiting hyperbolas H a ^r a \ about a are different for 
such a case. 

Our study has revealed the role of the input concentration Sin on the shape of the domain C. So we 
have computed numerically the best configurations (a*,r*) given by Proposition |4] as functions of Sin, as 
well as the corresponding output concentration S± (see Figure [7]) • The map Si„ M> a* given by (|28|l and 
(1301) being continuous, the discontinuity of the map a i-> f (a) leads to a discontinuity of the map Sin i V 5-^ . 
Consequently, there exists a threshold of 5j ra such that 

- below the threshold, the optimal buffered chemostat provides global stability, with performance close 
to the single chemostat i.e. 5* is close to A_ ; 

- above the threshold, the optimal stable configuration consists in a by-pass of the single chemostat 
without any buffer. The performance is significantly modified as 5* is larger than A+. 

According to Propositions [2] and this threshold corresponds to a value of Sin such that 5* (a*) = A + , 
where 5 is defined in (fT5|) . For values of 5j„ smaller than this threshold, the output concentration at steady 
state Si of the best configuration is thus bounded by the one computed for the limiting case when Si n get 
arbitrary close to the threshold (see Figure [8]). The values of 5j„ and 5* obtained at the threshold are given 
in Table [21 One can see on this example that the buffered chemostat allows a global stability for any value 
of Si n in the interval [0.777, 1.641] with an output at steady state less than 0.167, to be compared with the 
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limiting hyperbola when 
S(a)>A+ 




limiting hyperbola when 
S(oc)<A. + 

Figure 6: The limiting hyperbolas H a f ( a } about a = a (for Si n = 1.4). 



by-pass 



Figure 7: The best configuration (a*,r*) with S$, as functions of Si, 




fcSf X + Sin 

Figure 8: Determination of S$ for Si n at the threshold. 



Sin 


a* 


r* 




~ 1.641 


~ 0.543 


~ 0.561 


~ 0.167 



by-pass 



Table 2: Characteristics of the best configuration at the threshold value of Si n - 

value 0.103 of the locally stable equilibrium of the single chemostat (see also Figure[7]). 

In industrial applications, the attraction of the wash-out equilibrium is undesired because it presents a 
risk that may ruin the culture in case of disturbance, temporarily pump breakdown or presence of toxic 
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material that could drive the state in the attracting basin of the wash-out equilibrium. It imposes also to 
ensure that initial condition belongs to the attracting basin of the desired equilibrium. A common technique 
to overcome theses difficulties and allow an initial stage with a small concentration of biomass, is to control 
the input flow rate Q with a stabilizing feedback [51 |35] (it consists in finding a feedback law that reduces 
the flow rate when the state belongs to the attracting basin of the wash-out equilibrium) . But this solution 
requires an upstream storage and an actuator. The design of a buffered chemostat is thus an alternative 
that does not require any upstream storage nor feedback control. In real world applications, it may happen 
that the growth function /i(-) is not perfectly known or uncertain. Then choosing a buffered configuration 
not too close from the boundary of the domain C provides a robustness margin for the global stability. 

When the characteristics of the input flow cannot be changed, a simple solution consists in increasing 
the volume of the vessel, so that the dilution rate is small enough to ensure that condition of Case 3 of 
Proposition Q] is fulfilled. The relative increment AV/V has then to satisfy the condition 



S in <£{S>0\ (m(S) > 



1 



AV 
V 



that is equivalent to have 



AV 



1 



> 



1 



(34) 



(35) 



V ' (i(Si 7 , 

Note that under Assumptions Al and A2, this last number is positive. This solution increases significantly 
the residence time in the tank and induces additional financial costs. Instead of choosing a larger volume V, 
we show that adding a buffer can be an interested alternative to improve the stability of a given bioprocess. 
For the parameters given in Table [TJ we have compared numerically 

- the smallest relative increment of the volume of the single chemostat to be globally stable, given in 
13, 



- the smallest relative size of the buffer to be added for the buffered chemostat to be globally stable, 
given by Proposition [5] (that imposes to choose a = a*), 

as functions of the input concentration Si n (for values larger than A+ for which the bi-stability occurs with 
a dilution rate equal to one, cf Proposition [2]). One can clearly see on Figure [§] the advantage of the buffered 



AV 

V - single tank 





buffered 



single tank 

Sin 



Figure 9: Comparison of minimal increase of volume, and output nutrient concentration, to obtain the global 
stability. 



chemostat that requires less volume augmentation. The output concentrations are also drawn for both 
configurations with the minimal volume augmentation. According to Remark [2l these concentrations are 
always smaller than A_. This example demonstrates the flexibility of the buffered chemostat in the choice of 
possible configurations, with two parameters than can be tuned (see Figure [TO]) . while the single chemostat 
is penalized with only one parameter, requiring larger increments of volume and providing (too) low output 
concentrations. 
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V 1 ,~ S 1 5 1 .', 1 5 ' . ' . ' , ' , ' . ' . a " •• 1 1 1 1 1 a 1 1 1 1 a 1 1 u a " i 1 1 .:. 1 .:, 1 .;. 1 ,;, 1 ,;. 1 ,;, 1 .» 1 .:, a 

H(S in ) >l(Sto) ™ °" °" " "* °" "* °" °Jl(SiJ 

Figure 10: The sets of configurations (a, (3) that ensure global stability, to be compared with the minimal 
relative increase of volume of the single chemostat (for different values of Si„). 



Finally note that due to the robustness property that is obtained for the stability in the first tank when 
using a buffered chemostat, the presence of biomass at initial time is necessarily only in the buffer tank (see 
Proposition [3]) . This property possesses some advantages for the practitioners in industrial frameworks for 
the seeding phase. 

6 Conclusion 

The main message of the present study is that particular spatial structures can bring stability to unstable 
bioprocesses: 

- We have shown that a buffered interconnection of two volumes can globally stabilize the chemostat 
model when it exhibits a bi-stability, while preserving the same total volume and input flow, which is 
not possible with parallel or serial interconnections. 

- We have provided a characterization of the set C of all buffered configurations that enjoy this property, 
and study their performances at steady-state in terms of nutrient conversion. 

- Our study has emphasized the influence of the input concentration Si n on the shape of the set C and 
the design of the best configurations, that could exhibit a threshold on the value of Si n above it a 
by-pass is more efficient. More precisely, we have shown that this threshold can be computed from the 
function 

ip(s) = (S in - s)n(s) 
(where n(-) is the nutrient uptake rate assumed to be non-monotonic). 

- We have studied the minimal buffer volume to add to a single chemostat to obtain a global stability, and 
show how the flexibility of the buffered interconnection allows significantly less volume augmentation 
than increasing the size of the single chemostat. 

Those results provide new insights on the role of spatial structures in natural ecosystems, and new strategies 
for piloting bioprocesses designing volume and input flow of a buffer. Our study has considered single strain. 
According to the Competitive Exclusion Principle, it is not (generically) possible to observe more than 
one species at steady state in the buffer tank, but this does not prevent to have coexistence with another 
dominant species in the main tank, which is not possible with a single chemostat. Consequently, it might be 
relevant to study the performances of the buffered chemostat choosing a different strain in the buffer, that 
could be the matter of a future work. 

Appendix 

We recall a result from [551 Theorem 1.8] about asymptotic behavior of trajectories of asymptotically au- 
tonomous dynamical system. 
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Theorem. Let <I> be an asymptotically autonomous semi-flow with limit semi-flow 0, and let the orbit 
0$(t,£) have compact closure. Then the uj-limit set w$(t, £) is non-empty, compact, connected, invariant 
and chain-recurrent by the semi-flow 9 and attracts $(t, t, £) w/ien i — > oo. 
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